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Viscoelastic relaxation of topographic highs on Venus 
to produce coronae 
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Abstract. Coronae on Venus are believed to result from the gravitationally driven 
relaxation of topography that was originally raised by mantle diapirs. We examine 
this relaxation using a viscoelastic finite element code, and show that an initially 
plateau shaped load will evolve to the characteristic corona topography of central 
raised bowl, annular rim, and surrounding moat. Stresses induced by the relaxation 
are consistent with the development of concentric extensional fracturing common on 
the outer margins of corona moats. However, relaxation is not expected to produce 
the concentric faulting often observed on the annular rim. The relaxation timescale 
is shorter than the diapir cooling timescale, so loss of thermal support controls 
the rate at which topography is reduced. The final corona shape is supported by 
buoyancy and flexural stresses and will persist through geologic time. Development 
of lower, flatter central bowls and narrower and more pronounced annular rims and 
moats enhanced by thicker crusts, higher thermal gradients, and crustal thinning 
over the diapir. 


Introduction 

Ever since coronae were recognized as a unique land- 
form on Venus from images returned by the Venera 15 
and 16 spacecraft [Barsukov et ai, 1984, 1986; Niko - 
layeva et ai, 1986], their generally elevated topography, 
central bowl, tectonized raised annular rim, and sur- 
rounding moat have led many investigators to attribute 
their formation to endogenic processes, principally grav- 
itational modification of topography raised by mantle 
flow [Barsukov et ai, 1984; Pronin and Stofan . 1990; 
Stofan and Head, 1990]. This view was reinforced by 
the more extensive and detailed images and topographic 
data returned by the Magellan mission [Solomon et ai, 
1991; Squyres et ai, 1992a]. In this model, corona for- 
mation begins with mantle diapirism raising the surface 
to form a broad dome that may be radially fractured. 
Such domes were posited as a precursor to coronae [ Sio - 
fan and Head 1990], and several examples of radially 
fractured domes were revealed by Magellan [Squyres 
et ai, 1992a; Stofan et ai, 1992]. As the diapir flat- 
tens against the lithosphere, the model suggests that 
the dome assumes a more level, plateau shape [Janes et 
ai, 1992] which then gravitationally relaxes with diapir 
cooling to form the final coronal topography [Stofan et 
ai, 1991; Janes et ai, 1992]. 

Geophysical modeling has confirmed the ability of 
the proposed processes of diapiric uplift and subse- 
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quent relaxation to form the majority of topographic 
and tectonic characteristics of radially fractured domes 
and coronae [Stofan et ai, 1991; Janes et ai, 1992]. 
However, models to date have involved simplifying as- 
sumptions that have limited their usefulness in depict- 
ing the evolution of coronae. In particular, no model 
to date has addressed what the topographic end state 
of a corona should be. The studies of the relaxation 
portion of corona evolution by both Stofan et ai [1991] 
and Janes et ai [1992] employed the viscous model of 
Bindschadler and Parmeniter [1990]. This model in- 
cludes only viscous rheology and does not consider the 
elastic properties of the crust and mantle, which to- 
gether with isostacy, will ultimately support the load 
represented by the high-standing topography of coro- 
nae. The result of this simplification is that this model 
predicts that topography will ultimately relax entirely, 
producing aflat plain. In addition, although this model 
uses non-Newtonian flow laws for viscous rheology, it 
linearizes these laws by assuming that a constant, char- 
acteristic stress operates throughout the relaxation pro- 
cess, thus effectively treating viscosity as Newtonian. 
This approximation of viscous behavior permits pre- 
dicted stresses during relaxation to vary from a few tens 
of MPa to a few tens of thousands of MPa, depending 
on the exact choice of parameterizing constants [Stofan 
et ai, 1991; Janes et ai, 1992], and therefore fails to 
place reliable limits on the ability of the relaxation to 
produce the tectonism attributed to this stage of corona 
formation. 

Modeling of folding on Venus [Zuber, 1987; Zuber and 
Parmentier, 1990] indicates that the crust and mantle 
are both likely to have upper portions in which rheo- 
logical properties are dominated by brittle/elastic be- 
havior and lower portions in which rheology is dom- 
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inated by viscous behavior. The transition between 
dominant behaviors is governed by the brittle failure 
criterion and the power law creep behavior of the ma- 
terials involved, as well as the temperature and stress 
profiles with depth. Recent experimental work [Mack- 
well et ai, 1994] suggests that under the dry conditions 
prevalent on Venus, a diabase crust will have a much 
more viscous rheology than previously thought. The ef- 
fects of such a rheology have not yet been evaluated for 
their effect on viscous relaxation of topographic highs 
to produce coronae. 

In this paper, we address some of the shortcomings of 
previous modeling by employing a finite element analy- 
sis that includes elastic behavior, stress dependent non- 
Newtonian viscosity, and layered rheology. Our purpose 
is to answer several basic questions regarding the grav- 
itational relaxation phase of corona formation. These 


open questions include the timescale of relaxation, the 
final topographic state, and the stress state and pre- 
dicted fracture patterns at various times. We examine 
the dependence of each of these on crustal thickness, 
near-surface thermal gradient, choice of flow law for di- 
abase, and the initial geometry of the raised topography 
and that of the crust/mantle boundary. 

The purpose of this paper is to develop an under- 
standing of the nature of relaxation of topography on 
Venus to produce coronae and particularly to determine 
the final equilibrium topography and stress state. It is 
important therefore to keep in mind the general char- 
acteristics that any model of corona relaxation must be 
able to reproduce. Coronae vary greatly in size, but 
most have radii between 100 and 200 km [Stofan et ai, 
1992]. Mature coronae (Figure 1) are characterized by 
several topographic and tectonic signatures [Squyres et 



Figure 1 . Perspective view of Earhart, a corona located at 70°N, 137°E. View is from the SW. 
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Figure 2. Topographic profile across Earhart. Profile extends from NW to SE through the 
center of the corona. 


al, 1992a], In topographic profile (Figure 2; also see 
profiles of Stofan et al. [1992] and Janes et al. [1992]) 
coronae commonly show a central bowl-shaped or level 
area that stands a few hundred meters above the level 
of the surrounding plains. In a few examples, the lowest 
portion of the central area actually lies below the plains. 
The next feature outward is a raised annular rim 40 to 
100 km wide that lies 300-800 m above the plains and 
often shows evidence for concentric extensional fractur- 
ing. The rim falls off into a surrounding moat, typically 
30 to 50 km wide and with a floor several hundred me- 
ters lower than the surrounding plains. Many outer 
moat walls also also show concentric fractures, gener- 
ally attributed to extensional flexure and downwarping 
during moat formation [Stofan et al., 1991; Janes et 
al., 1992; Squyres et al., 1992a; Sandwell and Schubert, 
1992a]. We will demonstrate that most of these features 
will develop during the relaxation of plateau-shaped to- 
pographic highs on Venus, and that the topographic 
signature of central raised depression, annular rim, and 
encircling moat will persist through geologic time. 

Model 

TECTON 

To study the relaxation of topographic highs on Venus, 
we use the viscoelastic finite element code TECTON 
[Melosh and Raefsky, 1980]. We use an axisymmetric 
grid (Figure 3) in which the side nodes are fixed horizon- 
tally and are free to move vertically. The base nodes are 
fixed vertically and are free to move horizontally. All 
other nodes are free to move both horizontally and ver- 
tically. Circumferential movement is constrained by the 
axisymmetry to be zero. Since the base and and outer 
edge are fixed, they must be set far enough away from 
the area of interest that they do not affect its movement. 
Appropriate distances were found experimentally by in- 
creasing the depth and width of the grid until no change 
in results was observed with further increases. The grid 
used is 337 km deep and 1102 km wide, consisting of 
13,583 nodes and 14,269 elements. Each element is as- 
signed a set of rheological parameters including density, 


Young’s modulus, Poisson’s ratio, and non-Newtonian 
flow parameters as well as an initial stress state in three 
dimensions. In addition, forces proportional to any ver- 
tical movement are applied at certain nodes to simulate 
isostatic buoyancy. We take the density of the crust 
to be 3000 kg/m 3 and that of the mantle to be 3300 
kg/m 3 . Young’s modulus is taken to be 6 x 10 10 Pa, 
and Poisson’s ratio is 0.25 for both the crust and man- 
tle [ Turcotte and Schubert , 1982]. 

Failure Envelope 

TECTON is a viscoelastic finite element code, find- 
ing displacements and stresses by considering both the 
viscous and elastic rheological responses of materials. 
However, certain rheologic responses are either viscous 
or elastic, and the model starting conditions must re- 
flect these. Two such conditions are relevant to our 
model. One deals with the initial stress state through- 
out the model and the other with the application of 
buoyant restoring forces. To determine whether a given 
element is dominated by its elastic or its viscous re- 
sponse for these purposes, we adopt the sort of failure 
envelope, developed by Brace and Kohlstedi [1980], that 
has been used for Venus by a number of authors [e.g., 
Zuber , 1987; Banerdt and Golombek, 1988; Zuber and 
Parmentier t 1990]. 

Rocks can fail by either brittle failure (frictional slid- 
ing) or by plastic/ viscous failure (creep). Failure will 
occur by the weaker of the two possible mechanisms, 
generally by brittle failure above and by viscous flow 
below some critical depth. This critical depth occurs 
where the differential stress required for brittle failure 
equals the stress required to initiate viscous flow at a 
geologically significant rate. The critical depth at which 
the changeover between predominantly elastic (brittle) 
behavior and predominantly viscous behavior occurs is 
found by comparing the differential stress required to 
fracture rock as predicted by Byerlee’s law [ Byerltt , 
1978] with the differential stress required to cause the 
rock to flow at some geologically significant strain rate 
using experimentally determined flow laws. 

Byerlee’s law can be stated in compression as 
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Figure 3. (a) Finite element grid employed in this study, (b) Details of upper, on-axis portion 
of the grid. 


Acr c = 3.7(7, it. < 11 khar (1) 

Aa c = (2.1(7,) + 1.8 kbar a : > l 1 khar (2) 

where Acr c is the stress difference (i7j — <r 3 ) required to 
produce brittle failure in compression, and in extension 
as 

Act, = -0.79(7, it, < 5 3 khar (3) 

Act, = (-0.68(7,) + 0.57 khar a . > 5 3khar (-1) 

where Aa t is the stress difference required to produce 
brittle failure in extension. We take the vertical stress 
to be that due to the overburden 

k 

<r, = ^Tjp.gh,) (5) 

I = 1 

where g is gravity and p t and h t are the density and 
thickness of the overlying k rows of elements 
To model viscous behavior. non-Newtonian viscous 
flow laws are cast in the form: 

e = A exp(-Q//?r) ‘ (6) 

where i is the strain rate, c ^ is the differential stress, 
•ft is the gas constant, T is the absolute temperature, 
and A , n and Q are experimentally determined values. 
Equation (6) can be recast to determine the differential 


stress required for a given strain rate: 

(7 d = [(c/A) exp (Q/RT)) l ' n (7) 

We treat the temperature as following an error func- 
tion profile from a surface temperature of 735 K to a 
constant mantle temperature of 1700 K [Solomon and 
Head, 1991], The temperature profile is then fully de- 
scribed by defining the near-surface temperature gra- 
dient, or equivalently, the depth at which the mantle 
temperature becomes effectively constant with increas- 
ing depth. For purposes of determining the failure en- 
velope, we take the strain rate to be 1 0“ 15 s~ 1 . This is 
a geologically plausible strain rate, equivalent to a total 
strain of 3% over 1 m.y., and has been adopted in both 
terrestrial and venusian modeling [Bract and Kohlstedi, 
1980; Zuber, 1987; Banerdt and Golombtk , 1988; Zuber 
and Parmentier, 1990]. 

For both the crustal and mantle materials, above a 
critical depth at which the differential stress required 
to cause brittle fracturing and the differential stress re- 
quired to cause viscous flow are equal, the rheological 
response of the rocks will be dominated by elastic be- 
havior. Below this depth the rheology will be dominated 
by viscous behavior. Since we are concerned largely 
with the initial self-compression of the crust and man- 
tle, as well as with the downward flexure of the litho 
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Figure 4. Failure envelopes for the crust and mantle on Venus. The model shown is for an 11 
km thick diabase crust [ Cartstan , 1982] overlying an olivine mantle [ Goetze , 1978], a near-surface 
temperature gradient of 15 K/km, and a strain rate of 10“ 15 s _1 . Straight lines follow Byerlee’s 
law for brittle fracture, while curved lines follow the non-Newtonian flow laws. The weaker 
of the two mechanisms is indicated by the heavy solid lines. Brittle-ductile transition depths 
are indicated by dash-dot lines with elastic strength dominating above and viscous strength 
dominating below these lines. Dotted lines indicate viscous strengths for strain rates of 5 x 
1CT 15 s -1 and 2 x 10~ ^s' 1 , showing that the depth of changeover between dominant rheological 
behavior is relatively insensitive to the particular choice of strain rate. 


sphere under a load, we use the compressional elastic 
critical stress to determine the depth of changeover be- 
tween dominant behaviors. A typical failure envelope 
calculated in the fashion described above is shown m 
Figure 4. 

Initial Conditions 

Stress state. Changes with depth in the dominant 
rheological response have two important implications 
for the starting conditions of the finite element grid. 
The first of these is the initial stress state of the rocks. 
Rocks dominated by elastic behavior are assumed to be 
initially in a state of lithostatic stress 

0V = cr* = aj 3 (8) 

where c T is the radial stress and a * is the hoop stress, 
while rocks dominated by viscous behavior are assumed 
to be initially in a state of hydrostatic stress: 

(Tr = cr* = a 2 (9) 

We initially take the vertical stress, <7 i} to be that due 
to the overlying layers of material with no horizontal 
variation, i.e., we do not include the stresses due to the 
initial plateau load or any changes to the crust/mantle 
boundary below the load. To simulate these additional 
stresses, we make use of the initial elastic solution of 


TECTON, which returns any instantaneous displace- 
ments and resulting stresses due to the grid being ini- 
tially out of equilibrium. We take the resulting stresses 
for each element as the starting stress state for a new 
elastic solution. This iterative process is repeated until 
all initial elastic displacements are less than 1 m. 

Buoyancy Forces. The second important implica- 
tion in the change of the dominant rheological response 
mechanism is the imposition of buoyant restoring forces. 
It is sometimes erroneously thought that buoyancy re- 
quires a density contrast between materials. However, 
a solid material of p = 1.0 g/cm 3 immersed in water 
is subject to buoyancy. It is, in fact, neutrally buoyant 
specifically because of this force exactly counteracting 
the force of gravity. Buoyancy is a fluid response [Caih- 
les, 1975], so that this restoring force must be applied at 
each level at which a predominantly elastic material lies 
on top of a predominantly viscous material, or at the 
boundary between two immiscible fluids. In TECTON, 
such a force is applied at nodes and is equal to pAgbz 
where A is the area defined by the horizontal distance 
between the centers of the elements on either side of the 
node and (since TECTON treats axisymmetric grids as 
extending 1/2 radian from both sides of the plane of the 
grid for computational purposes) a concentric distance 
equivalent to 1 radian at the radial distance of the node, 
and bz is the total vertical displacement of the node. At 
the changeover from elastic to viscous crust the density 



21,178 


JANES AND SQUYRES: VISCOELASTIC RELAXATION OF CORONAE 


is that of the crust. At the changeover in the man- 
tle, or at the top of the mantle if the entire mantle is 
dominated by viscous rheology, the density is the differ- 
ence between the crustal and mantle densities, since the 
buoyancy of the crustal material will have already been 
applied at the appropriate level in the crust. For cases 
where the entire crust is dominated by elastic behavior, 
a single set of buoyancy forces is applied in the mantle, 
proportional to the mantle density. 

Time Steps and Ending Conditions 

The size of the time steps used during model runs is 
determined on the basis of the smallest Maxwell time 
for any element in the grid, defined as the ratio of that 
element’s viscosity to its Young s modulus. Time step 
size is constrained to be less than 10% of the Maxwell 
time during most of each run, but can be gradually in- 
creased to 30% of the Maxwell time at late stages when 
the relaxation is nearly complete and relaxation rates 
are small. Model runs using shorter time steps show no 
changes in calculated topography or stress. The limiting 
Maxwell time increases as viscous relaxation progresses. 
We have determined that the limiting Maxwell time as 
a function of time can be well fit by the relationship 

iog(r M ) = a + [6 logi] (10) 

where tm is the Maxwell time and t is the time. In prac- 
tice, we determine the initial limiting Maxwell time dur- 
ing the elastic iterations used to determine the starting 
stress state. We then use a time step size which is 10% 
of the initial Maxwell time and allow the grid to deform 
through a few thousand time steps. Actual Maxwell 
times are then determined at specific time steps, and 
the values of a and 6 in (10) are determined and used 
to calculate the predicted r m as a function of time for 
a particular run. Time step size is then increased as 
Maxwell time increases. Runs terminate when equilib- 
rium is reached and velocities have dropped to zero. At 
this point, remaining topography is supported by a com- 
bination of isostatic and elastic forces, and no further 
relaxation occurs. 

Crust/Mantle Boundary 

The process thought to produce coronae [Stofan et 
al., 1991; Janes et al. , 1992; Squyres et at., 1992a] 
begins with a rising mantle diapir which impinges on 
the underside of the lithosphere, raising domical topog- 
raphy which, as the diapir spreads beneath the litho- 
spheric lid, becomes more plateau shaped. The rising 
diapir also carries heat which would serve to raise the 
near-surface temperature gradient and to lower the crit- 
ical differential stress required for viscous flow in both 
the upper mantle and lower crust. If the crust is dy- 
namically coupled to the mantle flow at its base, lower 
portions of the crust may then become entrained in the 
outward mantle flow above the diapir head, creating 
a cavity in the crust which would be filled with rising 
diapir and/or mantle material [Bxndschadltr and Par - 
mentier, 1990]. Such crustal thinning would result in a 


higher density and Theologically separate mantle mate- 
rial in place of the lower portions of the thinned crust. 
We model such crustal thinning and infilling in some 
of our runs by assigning mantle rheology and density 
to elements under the load, extending upward from the 
mantle and laterally to a distance equal to the radius 
of the plateau. 4 

Results 

“Standard” Case 

Many studies have determined plausible ranges for 
Venus’ crustal thickness and thermal gradient. Deter- 
minations of crustal thickness from modeling of forma- 
tion of regularly spaced tectonic lineations [ Zuber , 1987; 
Banerdt and Golombek, 1988; Zuber and Parmentier, 
1990] yield values of 5 to 30 km. Crater relaxation 
models [Gnmm and Solomon , 1988] give 10 to 20 km, 
and inversion of gravity and topography [ Hemck and 
Phillips , 1992] implies a global average of 15 km. The 
crust could be significantly thinned over sites of man- 
tle upwelling, particularly if the mantle flow breaches 
the upper elastic mantle [ Btndschadler and Parmentier, 
1990]. To cover a meaningful range of values, we exam- 
ine crustal thicknesses of 5, 11, and 21 km. 

Thermal gradients have also been determined by many 
of the same studies [Zuber, 1987; Banerdt and Golombek, 
1988; Zuber and Parmentier , 1990; Gnmm and Solomon, 
1988] to range between 10 and 25 K/km. Simple con- 
duction models [Solomon and Head , 1982] and flexural 
studies [Solomon and Head , 1990] yield similar values. 
At the extremes, Sandwell and Schubert [1992b], mod- 
eling lithospheric flexure at large coronae, found ther- 
mal gradients as low as 7 K/km, while Squyres el al. 
[1992b], modeling buckling wavelengths in the plains, 
found values as high as 30 K/km. We examine near- 
surface temperature gradients of 5, 15, and 30 K/km. 

As a “standard” case for comparison with later re- 
sults, we. first examine the case of a 1 km high plateau 
load, 150 km in radius, on an 11 km thick crust with 
a near-surface temperature gradient of 15 K/km. We 
model the viscous behavior of the crust using the dia- 
base flow law of Caristan [1982] and of the mantle us- 
ing the olivine flow law of Goetze [1978]. For this case, 
the changeover from predominantly elastic to predomi- 
nantly viscous behavior occurs in the crust at a depth of 
680 m and in the mantle at a depth of 13 km (Figure 4). 
Relaxation predicted by the model (Figure 5a) proceeds 
very rapidly, with most of the motion occurring within 
the first year. The plateau relaxes to very near its final 
shape in approximately 100,000 years, with only lim- 
ited movement between 100,000 and 400,000 years. By 
about 400,000 years, equilibrium has been reached, the 
remaining load is supported by buoyancy and flexural 
stresses, and no further significant relaxation occurs. 

In reality, relaxation of topography raised by a hot 
diapir will be slower than our model predicts. This is 
because the model essentially assumes that the diapir 
that raised the topography cooled instantaneously, re- 
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Figure 5. (a) Topographic profiles through time for a 1-km plateau load on an 11-km crust 
with a near-surface temperature gradient of 15 K/km. Profiles are shown at 1 year, 10 years, 100 
years, etc. (b) Near surface stresses for the case shown in Figure 5a at 10,000 years. 


suiting in a state that is well out of equilibrium. This 
disequilibrium results in large stresses and tow viscosi- 
ties. In reality, relaxation rates would be controlled by 
the cooling rate of the underlying diapir rather than 
by the rheology of the mantle and crust. We can con- 
trast model-derived relaxation times with a character- 
istic cooling time for a diapir emplaced beneath the 
crust or lithosphere of Venus. Cooling of a diapir can 
be treated simply as a one-dimensional, time-dependent 
heat conduction case with a characteristic cooling time 
t c [Turcotte and Schubert, 1982] of: 

r c = / 2 A (11) 

where / is a characteristic length scale and k is the 
thermal diffusivity. As an illustrative case, consider 
a flattened diapir 35 km thick beneath a 5 km thick 
lithosphere, such as might be responsible for the raised 
plateau Selu [Janes et ai , 1992]. A length scale of 40 
km and a k of 1 x 10” 2 cm 2 s _1 [Turcotte and Schubert, 
1982] lead to a characteristic cooling time of 50 m.y. for 
the emplaced diapir, similar to the cooling timescale 
found by Stofan et ai [1991]. This cooling timescale 
is significantly longer, by 1 to 3 orders of magnitude, 
than any of the model-derived relaxation times that we 
obtain. Therefore the total time needed to completely 
relax to a final equilibrium state will be controlled by 
the cooling of the diapir rather than by any viscoelastic 
relaxation constraints. Similarly, relaxation will be ca- 


pable of responding essentially instantaneously to the 
cooling of the diapir. Any significant lag behind the 
withdrawal of thermal buoyancy will act to increase the 
differential stresses and thus decrease viscosity. This, in 
concert with viscosity decreases due to the higher tem- 
peratures in the vicinity of the cooling diapir, will keep 
relaxation rates high enough to follow cooling. We will 
therefore concentrate on the final topographic state of 
the relaxed corona, which is dependent on the rheology 
only, rather than the relaxation timescale, which is not. 

As expected, our model reproduces all the major to- 
pographic signatures of coronae. These form very early 
in the relaxation process and become more pronounced 
with time. The central bowl, high-standing rim, and 
exterior moat are all apparent early in the relaxation. 
Later, the central bowl sinks to lower elevations and 
becomes broader and flatter. The raised rim becomes 
narrower and more pronounced as its maximum eleva- 
tion decreases less rapidly than that of the central bowl. 
At its final equilibrium state, the crest of the rim sits 
approximately 300 m above the central bowl and 800 m 
above the surrounding plains level. Similarly, the moat 
becomes slightly narrower and deeper, reaching a final 
depth some 200 m below the plains. 

Stresses at a depth of 0.5 km (the midpoint depth 
of surface elements in the model) reach their maximum 
differential magnitude after approximately 10,000 years 
(Figure 5b), when the topography has relaxed to nearly 
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its final state. The stresses remain essentially the same 
with further relaxation. Stresses in the three princi- 
pal stress axes are compression al within the load, while 
the radial and hoop stresses become extensional beyond 
the load. The vertical stress is constant, representing 
pgh at a depth of 0.5 km within the upper layer of the 
load and the crust beyond the plateau load. Maximum 
differential stresses in the interior bowl region of the 
corona are compressional and relatively small, reaching 
a maximum of about 40 MPa (vertical minus radial) at 
a radial distance of 80 km. The vertical stress is the 
least compressional, and the radial stress is the most 
compressional, which could result in concentric thrust 
faulting if the stress levels were high enough [Ander- 
son, 1951]. However, rocks are strong enough in com- 
pression that the relatively low differential stress at this 
point would be unlikely to produce such faulting. In- 
deed, concentric thrust faults are not generally found 
within the central basins of coronae. 

Beyond the edge of the plateau, the near-surface ra- 
dial stress is the most extensional, while the vertical 
stress is the most compressional, leading to the expec- 
tation of concentric normal faulting. The maximum dif- 
ferential stress reached is 46 MPa at a radial distance of 
190 km. The differential stress is above 40 MPa over a 
range of 170 to 220 km from the corona center. This co- 
incides with the outer wall of the moat surrounding the 
raised rim, in the region where concentric grabens are 
observed on the moats of many coronae [S^t/yres et al., 
1992a]. The predicted stress level is close to the 50 MPa 
estimated to be required for fault propagation in Venu- 
sian surface rocks [ Janes et al., 1992]. There is a narrow 
area on the lowermost portions of the exterior wall of 
the moat where the radial stress is the most extensional 
and the hoop stress is the most compressional. Simple 
Anderson [1951] failure criteria would predict strike-slip 
faulting here if stress levels were high enough, though it 
is not clear that Anderson criteria can accurately pre- 
dict strike-slip faulting in such an environment [Schultz 
and Zuber, 1994]. No such faulting is generally observed 
here in coronae, and predicted differential stresses in 
this region are relatively small, approximately 30 MPa, 
so that in any case the absence of faulting is not unex- 
pected. The location of the maximum differential stress 
on the outer wall of the moat moves closer to the corona 
center with increasing relaxation. The maximum occurs 
initially at a radial distance of 230 km, moving inward 
as relaxation progresses. 

One striking discrepancy between model predictions 
and observed tectonics is the absence of a model pre- 
diction of concentric fracturing on the raised rim of the 
corona. Stresses on the rim are compressional, with the 
hoop stress the most compressional and the radial stress 
the least compressional so that, if anything, strike-slip 
faulting would be expected according to Anderson crite- 
ria. However, the differential stress in this region is only 
20 MPa or less. The lack of a prediction of the formation 
of concentric fractures on the rim during relaxation sug- 
gests that the fractures observed on corona rims form 
at some other stage of corona evolution For example, 


Koch [1994] suggests that this set of fractures forms 
during the transition from radially fractured domes to 
uplifted plateau, as bending of the surface is concen- 
trated at the edge of the plateau that will later form 
the raised rim. 

Crustal Thickness 

Proceeding from the “standard” case described above, 
we first examine the effects of crustal thickness on re- 
laxation. Figure 6a shows final profiles for crustal thick- 
nesses of 5, 11, and 21 km with a 15 K/km near-surface 
temperature gradient. It also shows a topographic pro- 
file across the corona Earhart. Relaxation takes longer 
for the 5-km crust than for an 11-km crust, reaching 
equilibrium after 800,000 years. This time span, how- 
ever, is still considerably shorter than the expected cool- 
ing time for the diapir which raised the initial plateau, 
so that relaxation time will still be controlled by cooling. 
The final topography is very similar to that produced 
by the thicker crust. 

For the 21-km-thick crust, all of the mantle is below 
the depth at which the dominant rheological behavior 
changes from elastic to viscous. The entire mantle is 
therefore initially in hydrostatic equilibrium, and man- 
tle buoyancy forces are applied at the top of the man- 
tle. The equilibrium topography is reached after only 
140,000 years of model relaxation. The topographic sig- 
natures are more pronounced, with a broader flattened 
region in the interior and narrower raised rim and sur- 
rounding moat. Late in the relaxation of this case, the 
center of the bowl rebounds slightly as load stresses 
transmitted by outer portions of the raised bowl dissi- 
pate with the continued relaxation of that portion of 
the corona. This rebound results in the center of the 
corona having a gentle dome rising a few tens of meters 
above the lowest level of the interior bowl. 

The model derived topography for all three crustal 
thicknesses generally fits that of the corona Earhart, 
in that they match the height of the central bowl and 
raised rim as well as the depth of the moat. However, 
the actual rim and moat are more separated and the 
slope between them shallower than that predicted by 
the model. This is almost certainly due to the fact that 
the model starts with an arbitrarily steep, step function 
plateau, whereas in reality, the slope of the plateau will 
likely be less pronounced. We will show, however, that 
relaxation of initial loads as broadly sloped as might be 
formed by thick extrusive flows will produce neither a 
central bowl nor raised rim. 

In Figure 6b, we plot the maximum radial-vertical 
differential stress, which would be responsible for con- 
centric tectonics, for each of the three different crustal 
thicknesses. The stress state reached with a 5-km-thick 
crust is very similar to that for the 11-km case, having 
a maximum of approximately 45 MPa on the outboard 
wall of the moat. Relaxation of a 21-km-thick crust 
leads to a more sharply localized rim and moat (Figure 
6a), and the maximum differential stress state is also 
more sharply defined. The differential stress reaches 
approximately 55 MPa at a distance of 160 km, closer 
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Figure 6. (a) Final equilibrium topographic profiles for a 1-km plateau load on 5- (dotted), 
II- (solid), and 21- (dashed) km crusts with a near surface temperature gradient of 15 K/km. 
The heavy dash-dot line is a topographic profile across the corona Earhart. (b) Maximum near 
surface stresses developed for the cases shown in Figure 6a. 


to the center than both the 5-km and 1 1-km crust cases. 
Differential stress exceeds 50 MPa over an annulus be- 
tween 160 and 170 km and exceeds 40 MPa between 
150 and 180 km. Therefore the width of the area of 
possible concentric fracturing narrows with increasing 
crustal thickness. Compressional stresses on the rim 
are also larger for thinner crusts, but still only reach a 
maximum of 50 MPa. 

Temperature Gradient 

To examine the effects of the near-surface thermal 
gradient on plateau relaxation, we ran additional mod- 
els with different thermal gradients. Figure 7a shows 
the results for an 11-km-thick crust with temperature 
gradients of 5, 15, and 30 K/km, as well as a topo- 
graphic profile across the corona Mama Alipa. With 
the shallower gradient, the changeover to viscous be- 
havior occurs deeper in both the crust and mantle. In 
the crust, the changeover depth is 760 m, while in the 
mantle it is 32 km. For the steep 30 K/km gradient 
the crustal changeover depth is 590 m, and the entire 
mantle is dominated by viscous behavior. 

Due to the temperature dependence of viscosity, re- 
laxation takes somewhat longer with the shallower tem- 
perature gradient, reaching equilibrium in 1 m.y., still 
shorter than most expected cooling times. The final 


topography is less pronounced than that found for the 
steeper 15 K/km gradient. The central elevation of the 
bowl is higher, and the rim and moat have shallower, 
broader slopes. The rim stands only 150 m above the 
central bowl, and the moat reaches a depth of only 180 
m below the surrounding plains. These differences are 
most likely due to the thicker portions of the crust and 
mantle that are dominated by elastic strength, and that 
therefore are able to flexurally support the load without 
failing by viscous flow. 

With a thermal gradient of 30 K/km, relaxation pro- 
ceeds very rapidly in the model and reaches equilibrium 
in only 10,000 years. Otherwise, the effects on final to- 
pography are generally similar to those produced by a 
thicker crust, with a narrowly defined rim and moat 
and a gentle rise in the central depression as relaxation 
reaches completion. 

As was the case with Earhart (Figure 6a), model de- 
rived rim height and moat depth provide a reasonable fit 
to that observed at Mama Alipa (Figure 7a), although 
a regional slope prevents a better fit to the flexure out- 
board of the moat. The topography of Mama Alipa 
shows a relatively sharp slope between the rim and 
moat, indicating that it might have been more sharply 
plateau shaped initially than was the case for Earhart. 
The depth of the central bowl implies a thermal gradi- 
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Figure 7. (a) Final equilibrium topographic profiles for a 1-km plateau load on an 11-km crust 
with a near surface temperature gradients of 5- (dotted), 15- (solid), and 30- (dashed) K/km. 
The heavy dash-dot line is a topographic profile across the corona Mama Alipa. (b) Maximum 
near surface stresses developed for the cases shown in Figure 7a. 


ent of at least 15 K/km, since shallower gradients would 
result in a shallower depression 

The effect of thermal gradient on the maximum dif- 
ferential stress developed during relaxation (Figure 7b) 
is similar to that produced by crustal thickness vari- 
ations (Figure 6b), but is more pronounced For a 
5 K/km gradient, the maximum differential stress is 
low and broadly distributed compared to the higher 15 
K/km case. In the region of concentric extension it 
reaches only about 25 MPa at a distance of 240 km. 
The central compressional stress is approximately 40 
MPa. Conversely, a steep 30 K/km gradient produces 
stresses that are both greater in magnitude and more 
narrowly confined to the moat and rim region Stresses 
reach a maximum differential stress of about 60 MPa in 
the area of concentric normal faulting on the outer wall 
of the moat at a distance of 160 km. 

In summary, both higher thermal gradients and thicker 
crusts produce more sharply defined corona topography 
and higher stresses. The center of the initial plateau 
typically subsides to approximately one half its original 
height of 1 km. The exception is when a very low ther- 
mal gradient results in such thick elastically dominated 
portions of the crust and mantle that the center remains 
high. In the most advantageous case of a thick crust 
and a high thermal gradient, the differential stresses on 
the outer wall of the moat are greater than the 45 MPa 


maximum brittle extensional strength of the crust (Fig- 
ure 4) and the 50 MPa fault propagation strength of the 
Venusian crust inferred by Janes et al [1992]. 

Viscous Flow Law 

Recent laboratory experiments have indicated that 
under the very dry conditions on Venus, diabase may 
exhibit a more viscous behavior than previously thought 
[Mackwell et ai, 1994]. We have run models using this 
new flow law to determine the effect it would have on 
relaxation of plateaus on Venus. Figure 8a shows re- 
sults for 5-, 11-, and 21-km thick crusts with a viscous 
rheology that follows the flow law of Mackwell et al. 
[1994] and a near surface temperature gradient of 15 
K/km. For this rheology and temperature gradient, the 
entire 5 km thick crust is dominated by elastic behav- 
ior, while the 11- and 21-km thick crusts change from 
dominantly elastic to dominantly viscous at a depth of 
7 km. In all cases, the mantle changeover occurs at a 
depth of 13 km. Topographic profiles after relaxation 
are essentially unchanged from those obtained using the 
flow law of Caristan for the 5- and 11-km-thick crusts, 
implying that, at least when the crust is relatively thin, 
relaxation is controlled primarily by the mantle viscos- 
ity rather than that of the crust. However, there is a 
significant change in the final relaxed profile for thicker 
crusts. The 21 km thick crust using the Mackwell flow 
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Figure 8. (a) Final equilibrium topographic profiles for a 1-km plateau load on 5- (dotted), 11- 
( solid), and 21- (dashed) km crusts with a near surface temperature gradient of 15 K/km and 
employing the viscous flow law of Mackutll et al. [1994]. (b) Maximum near surface stresses 
developed for the cases shown m Figure 8a. 


law does not develop the low central uplift s#*en with 
the softer material. In addition, tlu* rim and moat are 
less sharply peaked, resulting in a final profile which is 
virtually indistinguishable from the thinner crust cases 

The more viscous crustal rheology similarly affects 
surface stresses (Figure 8b). Stresses for the .> and 
11-km-thick crusts are similar to those obtained using 
the Caristan flow law- (Figure 6a) However, stresses 
developed in a 21-km-thick crust are marked!) changed 
when using the more viscous Mack well flow law As was 
the case with topography, stresses remain essentially 
unchanged as the crust thickens and do not develop 
sharply peaked differential stresses in the area of the 
rim and moat. 

The main effect of using the more viscous diabase flow- 
law is to increase the time required for the models to 
reach equilibrium. However, even these extended limes 
are shorter than the expected diapir cooling times. 
Maximum differential stresses become slightly larger, 
develop over a wider area, and occur significantly later 
in the relaxation process. For crusts that are as thick 
as 21 km, the increased crustal viscosity also produces 
a broadening of the moat and nm. 

Crust /Mantle Boundary 

The models we have examined so far have assumed 
that the crust/mantle boundary occurs at the same 


depth under the raised topography as under areas be- 
yond. However, the generalized model of corona forma- 
tion posits that the raised topography is due to uplift of 
the elastic lithosphere by a rising mantle diapir. Mantle 
flow above the diapir will be outwards along the base 
of the crust, raising the possibility that the lower crust 
will become entrained in this flow, resulting in crustal 
thinning [Btndschadler and Parmentier, 1990]. In or- 
der to examine the effects of replacing these portions of 
the crust with mantle material, we have run models in 
which elements that lie at depths between the nominal 
crust/mantle boundary and a depth of approximately 
one half the crustal thickness, and which are under the 
initial plateau topography, are assigned mantle material 
properties. 

The replacement of crust with mantle material results 
in several changes in final topography, and maximum 
differential stress. For a 5 km crust (using the Caristan 
flow law) with mantle material replacing the lower 2 km 
of the crust, the excess mass causes the central bowl to 
descend an additional 200 m to a height 400 m above 
the plains level (Figure 9a). The rim and moat topog- 
raphy becomes narrower and more pronounced, the rim 
standing 350 m above the bowl and 750 m above the 
plains, while the moat reaches a depth of nearly 250 
m. For an 11-km-thick crust with mantle material fill- 
ing the lower 5 km below the plateau, the bowl of the 
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Figure 9. (a) Final equilibrium topographic profiles for a I-km plateau load on a 5-km crust 
with mantle material intruding into the lower 2 km of the crust under the load (dotted), an 11-km 
crust with mantle material intruding into the lower 5 km (solid), and a 21-km crust with mantle 
material intruding into the lower 10 km (dashed). All cases have a near surface temperature 
gradient of 15 K/km (1>) Maximum near-surface stresses for the cases shown in Figure 9a. 


corona would now lie some 300 m above the surround- 
ing plains. The rim lies 400 m above the center and 700 
m above the plains, while the moat is approximately 
300 m deep. The rim and moat topography is also sim- 
ilar to the 5 km case, being narrower and steeper than 
was the case without crustal thinning (Figure 5a). In 
the case of a 21-km-thick crust with mantle material 
infiltrating the lower 10 km. the central howl descends 
to nearly the level of the surrounding plains, lying only 
100 m above them (Figure 9a) This greater depression 
of the central bowl makes the rim more pronounced It 
now stands some 500 m above the center even though 
it lies only 600 m above the level of the plains. The 
moat becomes both deeper and wider, reaching a depth 
of 400 m and a width of about 50 km The rebound 
of the center to form a gentle dome which was evident 
in the 21-km-thick crust (Figure 6a) still occurs but is 
much less pronounced. 

Stresses, as well as topography, become higher and 
more sharply focused with replacement of the lower 
crust by mantle material. The peak differentia] stress on 
the outer wall of the moat for the 5 km thick crust is 53 
MPa (Figure 9b). Differential stresses are greater than 
50 MPa over radial distances of 175 to 200 km from 
the corona center. Without mantle material displac- 


ing portions of the crust, stresses in this region never 
exceed 50 MPa. Thus the replacement of crustal ma- 
terial with higher density mantle makes development 
of a broad annulus of concentric extensional fracturing 
more likely than would otherwise be the case. Com- 
pressional stresses within the bowl are also higher but 
reach a maximum differential stress of only 47 MPa, 
probably inadequate to produce compressional tecton- 
ism. Differential stresses in the 11 km thick crust reach 
a maximum of 64 MPa at a radial distance of 190 km. 
Stresses are greater than 50 MPa from 160 to 220 km, 
while compressional stresses in the bowl reach a maxi- 
mum of 60 MPa at a distance of 85 km. For the 21 km 
thick crust having the lower 10 km replaced with mantle 
material, stresses on the outer wall of the moat reach a 
maximum of 92 MPa and exceed 50 MPa from 160 to 
225 km. In this case, however, compressional stresses 
within the bowl are also much larger, reaching 71 MPa 
at a radial distance of 79 km, large enough that one 
might begin to expect to see concentric thrust faulting 
and/or folding under these conditions. 

Topography and tectonics similar to those our mod- 
eling predicts to result from significant crustal thin- 
ning are seen at a few coronae on Venus. For exam- 
ple, Squyres ei ai [1992a], describe the rim of the 
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Figure 10. Final equilibrium topographic profiles for a 1-km parabolic load on an 11-km crust 
with a near surface temperature gradient of 15 K/km. 


corona Ba’het as lying less than 400 m above the sur- 
rounding plains while its interior is at or slightly be- 
low plains level. They also note that, while the ma- 
jority of the concentric tectonics around Ba’het are ex- 
tensional in nature, the innermost ring of features at 
its northern end are compression al. Similarly, at least 
two other coronae exhibit compressional rings that are 
always the innermost of the circumferential tectonics, 
Ekajata [Squyres et al, 1992] and Sith [Stofan et al , 
1992]. 

Load Shape 

Previous viscous modeling of topographic relaxation 
on Venus [Stofan et al., 1991] indicated that an initially 
uncompensated plateau shape was required in order to 
form the characteristic topography of coronae. We have 
implicitly assumed that this is the case in the modeling 
shown so far. To test this assumption, we have run a 
model in which the initial topography has a cross sec- 
tion that would result from an outflow of a viscous fluid 
onto the surface. To determine the initial load shape, 
we adopt the similarity solution of Huppert [1982], who 
found the height of an axisymmetric outflow, as a func- 
tion of radial position and time, to be given by 

V .«) = Sr 3 (3 Qv/g') 1 '* < ( °- ,)/4 <HUZr) (12) 

where h is the height, r is the radial distance from the 
flow center, R is the radius of the flow, t is time, Q is the 
constant of proportionality between the volume of the 
flow and f°, v is the viscosity, g ' is the reduced gravity 
(equal to gAp/p), and a is a nonnegative constant. At 
times longer than the emplacement time of the flow 
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solving for these constants from the boundary condi- 
tions at the edge of the flow, 


Sr = Q (19) 

and at the center of the flow, 

/»(r.=0) = SR 3 CtWS(r= 0 )/U) (20) 

Figure 10 shows the topography of this kind of load, 
initially 1 km in height, during relaxation. At no time 
during the process does the characteristic profile of a 
corona, with a central depression and an annular rim, 
form. Only the exterior moat develops as the litho- 
sphere flexes under the imposed load. The inference 
drawn by Stofan et al. [1991] that the initial topography 
must be flat-topped and steep-sided therefore is robust 
to the addition of elastic rheology to relaxation model- 
ing. Such plateau-shaped topography should result as 
the diapir that produced the initial radially fractured 
dome flattens against the underside of the lithosphere 
[Janes et al. , 1992; Koch , 1994], and is also likely to 
be the source of concentric extensions! fracturing con- 
centrated at the edge of the plateau [Koch, 1994] which 
will later form the raised annulus. 

Discussion and Conclusions 

Finite element modeling of gravitationally driven vis- 
coelastic relaxation on Venus indicates that the topog- 
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raphy associated with coronae can develop over rela- 
tively short time spans and that actual relaxation rates 
will therefore be controlled by the cooling rate of the 
underlying diapir. This topography is then supported 
by a combination of buoyancy and flexural stresses and 
will persist through geologic time. Stresses on the outer 
wall of the flexural moat are in the proper orientation 
and, in some cases, of sufficient magnitude to produce 
the concentric extensional fracturing that has been ob- 
served. However, the stresses on the annular rim will 
not produce the fracturing that is generally observed 
there. A thicker crust and/or a higher near-surface tem- 
perature gradient will produce a narrower raised rim 
and moat, and a flatter interior bowl. Thinning of the 
crust and infilling with mantle material produce a lower 
interior bowl and deeper moat, while making the raised 
rim broader and higher relative to the depressed inte- 
rior. All three conditions, thick crust, high temperature 
gradient, and localized crustal thinning due to mantle 
flow, produce higher differential stresses on the outer 
wall of the moat, increasing the probability of the for- 
mation of concentric extensional fracturing. The choice 
of specific non-Newtonian flow law does not significantly 
affect corona formation, since relaxation will still be 
rapid enough that it can instantaneously respond to 
the loss of buoyancy as the diapir cools. A necessary 
starting condition to the production of coronae is a flat- 
topped plateau. The profiles that develop from our ini- 
tial plateau shape show more sharply peaked rims than 
are observed. This difference is most likely due to our 
initial plateau being very sharply defined, essentially 
as a step function, rather than having more rounded 
shoulders with shallower slopes 
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